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摘要 : 割 胶 促 进 橡胶 树 合成 天 然 橡胶 与 激活 乳 管 细胞 的 薪 莉 酸 信号 途径 密切 相关 。 茉 莉 酸 信号 途径 关键 环 


节 基 因 的 表达 水 平和 橡胶 生物 合成 相关 基因 的 表达 水 平 与 干 胶 产 量 正 相 关 。 但 这 些 基因 彼此 间 的 表达 相关 


性 尚 不 清楚 。 本 文通 过 qPCR 获得 


HbJAZ2、HbJAZ3、HbMYC1、HbMYC2、HbMYC3、HbMYC4、HbMYC5 和 6 个 橡胶 生物 合成 相关 


了 正常 割 胶 条 件 下 , 9 个 茉莉 酸 信号 途径 关键 环节 基因 HbCOIH « HbJAZI 
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大 | HbHRT2. 


HbSRPP. HbREF. HbHMGRI1. HbHRTI. HbGAPDH 在 5 个 橡胶 树 魏 克 汉 种 质 和 5 个 1981’IRRDB 种 质 胶 
乳 中 的 表达 数据 ; 通过 皮尔 逊 相关 系数 分 析 了 这 15 个 基因 彼此 间 的 表达 相关 性 , 分 别 获得 105 对 基因 的 双 
变量 相关 系数 (120 和 偏 相 关系 数 (r12.3)， 所 有 105 对 基因 的 双 变 量 相关 系数 C | r12 | ) 和 偏 相 关系 数 


一 


供 理论 基础 。 


据 ， 可 用 于 挖掘 、 筛 选 和 预测 与 橡胶 4 


关键 词 : 巴西 橡胶 树 ， 割 胶 ， 茉 
中 图 分 类 号 : Q943.2 


rl2.3 | ) 的 平均 值 分 别 为 0.486+0.220 和 0.304 € 0.211, 达到 0.05 的 差异 显著 水 平 。 其 中 , r12 与 T12.3 
方向 相同 的 63 对 (60%)， 方 向 相反 的 42 对 〈409% ); 
r12.3 | 的 82 对 (78.095%); 双 变 量 相关 显著 虱 


| r12 | < |rl2.3 | 的 23 X} (21.905%), | rl2 | > 


E P-0.05 的 76 对 (72.38%), P<0.01 If] 59 对 (56.19%); 
偏 相关 显著 性 P<0.05 的 21 对 (20%), P<0.01 的 16 对 (15.24%)。 结 果 显 示 ， 这 两 条 途径 中 的 基因 表达 
彼此 相关 ， 这 为 基因 表达 谱 分 析 中 普遍 采用 的 前 提 假 设 “ 功 能 相关 基因 的 表达 相关 ”提供 了 进一步 的 实验 证 
居 E 物 合成 相关 的 未 知 基因 ， 为 研究 橡胶 树 产 量 形成 的 分 子 调控 机 理 提 


莉 酸 信号 ， 橡 胶 生 物 合成 ， 基 因 表 达 ，Pearson 相关 
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Abstract: Tapping-enhanced rubber biosynthesis is closely related to the activation of jasmonic acid signaling in 
laticifer cells of rubber tree. The expression level of genes related to jasmonic acid signaling pathway and rubber 
biosynthesis were both positively correlated with the dry rubber yield. However, the exact relationship between the 
expression of gengs involved in jasmonic acid signaling and rubber biosynthesis is yet to be elucidated. In the 
present study, qPCR was used to analyze the expression of 9 jasmonic acid signaling genes, HbCOII, HbJAZI, 
HbJAZ2, HbJAZ3, HbMYCI, HbMYC2, HbMYC3, HbMYC4, HbMYCS5 and 6 rubber biosynthesis genes, HbHRT2, 
HbSRPP, HbREF, HbHMGRI, HbHRTI, HbGAPDH, in laticifer cells of 5 Wichham germplasms and 5 
1981’IRRDB germplasms following tapping them with S/2D d/3 tapping system, Hb18S was used as internal 
reference gene. The correlation between these 15 genes was analyzed by Pearson correlation coefficient. The 
bivariate correlation coefficient (r12) and partial correlation coefficient (112.3) of 105 gene pairs were obtained 
with the mean + standard deviation were 0.486 + 0.220 and 0.304 x 0.211, respectively. Among them, 63 gene 
pairs (6096) with r12 in the same direction as r12.3 and 42 genes pairs (40%) in the opposite direction; In the 
degree of correlation, 23 gene pairs (21.905960) with | r12 | less than | r12.3 |, and 82 gene pairs (78.095906) with | 
r12 | greater than | r12.3 |. There 76 gene pairs (72.3896) bivariate correlation coefficient (r12) were significant 
correlation at P «0.05 level and 59 gene pairs (56.19%) at P —0.01 level. In contrast, there less partial correlation 
coefficient (r12.3) were significant correlation at P<0.05 level (21pairs, 2096) and P<0.01 level (16 pairs, 
15.2496). These results suggest that the expression of genes involved in these two pathways were related to each 
other. And provides a further evidence for the assumption that “expression correlation of functional related genes", 
which is widely adopted in gene expression profile analysis. It can be used for excavating, screening and 
predicting unknown genes related to rubber biosynthesis, as well as provides a theoretical basis for studying the 
molecular regulation mechanism of rubber yield formation in Hevea brasiliensis Muell. Arg. 

Key words: Hevea brasiliensis, tapping, jasmonate signaling, rubber biosynthesis, gene expression, pearson 
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茉莉 酸 是 植物 逆境 响应 的 关键 激素 信号 (Qi etal., 2011) , COIN, JAZ 和 MYC 是 茉莉 
酸 信号 途径 的 关键 环节 (Chini et aL, 2009); JA-Ile, COIL, JAZ 和 MYC 通过 精细 的 互 作 
CDonnell et al., 1996; Chini et al., 2007; Fonseca et al., 2009; Qi et al., 2011 ) 参 与 植物 次 生 代 谢 
的 调节 (Qi et al., 2011; Schweizer et al., 2013) 。 巴 西 橡胶 树 中 的 天 然 橡胶 生物 合成 是 一 种 典 
型 的 植物 类 异 戊 二 烯 代谢 ， 受 茉莉 酸 信号 途径 的 调节 (Deng et al., 2018) ; 橡胶 生物 合成 途 
径 包 括 一 系列 橡胶 生物 合成 关键 酶 (Koyama & Tanaka., 1996) ， 如 REF (Dennis & Light, 
1989) ~ HMGR (Chye et al., 1992) ~ HRT CAsawatreratanakul et al., 2003) 、SRPP (Collins, 
2009) 等 ;茉莉 酸 信号 途径 对 这 些 酶 的 精细 调节 机 制 尚 不 完全 清楚 。 橡 胶 树 种 质 PR107、 
RRIM600、 热 明 628. WE 525 和 热 明 523 的 干 胶 产 量 普遍 高 于 RO/CM/10 44/160. MT/IT/13 
29/8. RO/C/8 24/104, RO/I/103 107 和 RO/CM/10 44/454， 并 且 莱 莉 酸 信号 途径 关键 环节 基 
HbCOII. HbJAZI. HbJAZ2、 HbJAZ3. HbMYCI. HbMYC2. HbMYC3. HbMYC4、 HbMYCS 
的 表达 和 橡胶 生物 合成 酶 基因 HbHRT2 . HbSRPP . HREF, HbHHMGRI , HbHRTI , HbGAPDH 
的 表达 与 他 们 的 橡胶 产量 正 相 关 《〈 杨 署 光 等 ，2019a; 杨 署 光 等 ，2019b) 。 但 是 ， 这 些 基 因 
表达 的 株 间 差 异 以 及 这 些 基 因 间 的 表达 相关 性 尚 不 清楚 。 本 研究 通过 分 析 这 些 基 因 在 这 10 
份 种 质 中 的 表达 相关 性 ， 可 为 研究 橡胶 树 产 量 形成 的 分 子 调控 机 理 提供 理论 基础 。 


1 材料 与 方法 
1.1 材料 


| 


10 份 巴西 橡胶 树 CHevea brasiliensis) 种 质 、 试 剂 、 耗 材 与 先前 的 报道 相同 〈 杨 署 光 等 ， 
2019a; 杨 署 光 等 ，2019b) 。 其 中 ，RO/C/8 24/104 为 2 Yk. RO/V103 107 为 4 株 ， 其 他 种 
质 各 3 株 ， 共 30 株 树 。 

12 方法 

1.2.1 材料 处 理 
生产 中 按 S/2D d3 正常 制 胺 ， 在 8 月 份 的 某 次 割 胶 时 ， 分 别 收集 前 10 min 流出 的 胶乳 ， 

于 提取 胶乳 总 RNA。 

1.2.2 总 RNA 的 提取 与 cDNA 的 合成 

用 曾 日 中 等 的 方法 《〈 曾 日 中 等 ，2003) 提取 胶乳 总 RNA。cDNA 第 一 链 的 合成 与 先前 

的 报道 相同 ( 杨 署 光 等 ，2019a; 杨 署 光 等 ，2019b) 。 

1.2.3 基因 表达 分 析 

1.2.3.1 qPCR 引物 合成 

引物 为 序列 来 自 文献 (TIAN et ab, 2013; 杨 署 光 等 ，2019a; 杨 署 光 等 ，2019b) ， 委 
托 Invitrogen 公司 合成 。 
1.2.3.2 qPCR 反应 

qPCR 反应 按 先 前 报道 的 方法 完成 〈 杨 署 光 等 ，2019a; 杨 署 光 等 ，2019b) 。 

1.3 数据 处 理 

以 Hb18S 作 内 参 ， 按 先前 报道 的 方法 〈 杨 署 光 等 ，2019a; 杨 署 光 等 ，2019b) 计算 目 

的 基因 的 相对 表达 值 。 
用 SPSS 软件 分 析 基 因 表 达 值 的 Pearson 相关 性 ，P<0.01 表示 极 显著 相关 ，P 二 0.05 
表示 显著 相关 ; 相关 性 的 密切 程度 参照 文献 报道 ( 朱 婉 丽 ,， 2019) 的 标准 , 根据 相关 系数 OD 
的 数值 大 小 ， 分 为 四 个 等 级 : 0< | r | 0.3 的 微弱 相关 ; 0.3< | r | 0.5 的 低 度 相关 ;0.5 
< |r| «0.8 的 中 度 相 关 ; 08 Er] <1 的 高 度 相关 。 

分 别 用 SPSS 软件 的 Duncan 检验 和 Excel TTEST (Array 1，Array2，Tails 2, Type 2) 
完成 多 重 比 较 分 析 和 双 样 本 比较 分 析 : P<0.01 表示 组 间 差 异 极 显 著 ，P<0.05 表示 组 间 差 
异 显 著 。 

用 GraphPad Prism 5 作 图 。 
2 结果 与 分 析 
2.1 橡胶 树种 质 胶乳 中 橡胶 生物 合成 调控 相关 基因 的 表达 分 析 

qPCR 分 析 表 明 ， 不 同 种 质 中 的 基因 表达 差异 明显 ， 相 同 种 质 的 生物 学 重复 间 也 能 达到 
极 显 著 差 异 ， 并 且 高 产 种 质 中 的 基因 表达 普遍 高 于 低产 种 质 (图 1、 图 2)。 在 高 低产 种 质 组 
iR, COH (图 1: AX. JAZI (图 1: B), JAZ2 (图 1: C), JAZ3 (图 1: D), MYCI (图 1: 
E), MYC2 (图 1: F)、MYC3 (图 1: G)、MYC5 (图 2: D. GAPDH (图 2: J), HMGRI 
(图 2: K), REF (图 2: MD. HRT2 (图 2: O) 差异 极 显 著 ， HRTI1 (图 2: ND 差异 显著 ， 
SRPP (图 2: L) 差异 不 大 ，MYC4 (图 1: H) 差异 不 显著 ; 前 后 者 之 比分 别 为 1.74、1.84、 
2.17、3.73、2.58、2.38、8.72、0.34、2.58、1.36、2.28、3.02、1.84、1.36、1.01。 总 的 来 说 ， 
15 个 基因 在 高 低产 种 质 组 间 的 平均 表达 差异 极 显 车， 两 者 之 比 为 2.08 (图 2: PO. 
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E: 横 坐 标 ， 从 左 到 右 〈1-10) 依次 代表 魏 克 汉 种 质 PRIO7. RRIM600. 4&5 628. HE 525. ME 523 以 
及 1981'IRRDB 种 质 RO/CM/10 44/160. MT/IT/13 29/8. RO/C/8 24/104, RO/I/103 107. RO/CM/10 44/454; 
不 同 大 写字 母 (或 **) 表示 组 间 差 异 极 显著 (P<0.01), 不 同 小 写字 母 (或 *) 表示 组 间 差 异 显著 (P<0.05)。 
每 根 柱子 代表 一 个 生物 学 重复 。 下 同 。 
Note: X-axis, from the left to right (1 to 10), the sample are Wickham germplasms PR107, RRIM600, Re ken 628, 
Re ken 525, Re ken 523, and 1981 'IRRDB germplasms RO/CM/10 44/160, MT/IT/13 29/8, RO/C/8 24/104, 
RO/T/103 107, RO/CM/10 44/454; The different capital letters (or **) indicate significance at P 0.01, the 
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different small letters (or *) indicate significance at P<0.05. Each column represents a biological repetition. The 
same below. 
1 橡胶 树种 质 胶 乳 中 橡胶 生物 合成 调控 相关 基因 的 表达 分 析 


Fig.1 Expression analysis of genes related to rubber biosynthesis regulation in latex from rubber 


tree germplasms 
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: 每 个 生物 学 重复 的 相对 表达 量 为 15 个 基因 的 平均 数 + 标 准 差 (P: AVERAGED. 
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Fig.2 Expression analysis of genes related to rubber biosynthesis regulation in latex from rubber 


Note: Data are x + s of 15 genes in each biological repetition (P: AVERAGE). 


图 2 橡胶 树种 质 胶乳 中 橡胶 生物 合成 调控 相关 基因 的 表达 分 析 


tree germplasms 


2.2 橡胶 树种 质 胶乳 中 橡胶 生物 合成 调控 相关 基因 表达 的 相关 性 分 析 


本 研究 共 分 析 了 15 个 


因 间 的 表达 相关 性 ， 分 别 获 得 105 对 基因 的 双 变 量 相关 系数 


(r12) 和 偏 相关 系数 (0012.30 (图 3、 图 4)。 其 中 ，r12 5 r12. 3 方向 相同 的 占 60%， 相 
反 的 占 40%; 微弱 相关 、 低 度 相 关 、 中 度 相 关 、 高 度 相 关 的 r12 分 别 为 20.952%、 33.333%、 


37.143%、8.571%， 在 r12. 3 中 分 别 为 54.286%， 
< | 112.3 | 的 为 21.905%， 
为 72.38% 和 56.1996, 1E r12. 3 H 
途径 相关 基因 表达 的 平均 值 与 6 个 橡胶 9 
极 显著 正 相 关 。 
分 析 了 单个 基 


105 对 基 


因 分 别 与 其 余 14 个 基因 


22.667%, 18.095906, 0.952%; | rl2 | 
| r12 | > | rl2.3 | 为 78.095%; 显著 相关 和 极 显著 的 r12 分 别 
FP 分 别 为 20% 和 15.2496; Æ 30 个 植株 中 ，9 个 茉莉 酸 信 和 号 


FE 物 合成 相关 基因 表达 的 平均 值 间 的 r12 为 0.892, 


间 Pearson 相关 系数 的 平均 值 (图 4: P)。 所 有 
和 | r12.3 | 的 平均 值 分 别 为 0.486 和 0.304, 两 者 差异 显著 (图 AP: Total); 


HRT2、REF、COI1、JAZ1、JAZ3、MYC1、GAPDH、MYC2、MYC3 的 1r12| 与 1r123| 差 
异 显 著 ， 而 HMGRI. SRPP. HRTI. JAZ2. MYC4. MYCS 的 | rl2 | 5 |1123 | 差异 不 显著 ; 


各 组 中 | rl2 | 普遍 大 于 1023 | ， 但 彼此 间 的 差异 显著 性 
fi, MYCL2,3 与 其 他 14 个 基 
员 中 ，JAZ1,2,3 与 
JF; Xs. REF. GAPDH, MYCI. MYC2 fi 
MYC4 和 MYC5. 


L 有 广泛 的 交集 ;5 个 MYC 家 族 成 


因 的 相关 性 相当 ， 


的 相关 性 差异 不 大 


极 显 著 高 于 MYC4,5; 3 个 JAZ 家 族 成 
， 仅 JAZL2 间 的 |rl2 | 达到 0.05 的 显著 
Il MYC3 与 其 他 14 个 基因 的 相关 性 高 于 
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Note: ** indicates significant correlation at P<0.01, * indicates significant correlation at P=0.05(2- 
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Significance). The same below. 


图 3 橡胶 生物 合成 调控 相关 基 


Fig.3 Pearson correlation coefficient among the expression of genes related to rubber biosynthesis 
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Fig.4 Pearson correlation coefficient among the expression of genes related to rubber biosynthesis 
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3 讨论 与 结论 


橡胶 树 的 胶 生物 合成 是 一 种 典型 的 植物 类 异 戊 二 烯 代 谢 ， 茉 莉 酸 信号 途径 与 类 异 成 二 
烯 代谢 途径 密切 相关 (Deng et al., 2018) ， 橡 胶 树 中 这 两 条 途径 的 相关 基因 表达 与 橡胶 产量 
正 相 关 《〈 杨 署 光 等 ，2019a; 杨 署 光 等 ，2019b; 卢 世 香 ，2010; Æ, 20130 。 该 研究 通过 
增加 生物 学 重复 ， 进 一 步 证 实 了 这 些 基 因 的 表达 与 橡胶 产量 的 相关 性 ; 同时 ， 该 研究 结果 显 
示 ， 这 两 条 途径 中 的 基因 表达 彼此 间 有 具有 一 定 的 相关 性 ,说 明 功 能 相关 的 蛋白 ,基因 表达 倾 
向 相关 ; 这 为 基因 表达 谱 分 析 中 普遍 采用 的 前 提 假 设 “ 功 能 相关 基因 的 表达 相关 ”提供 了 理论 
依据 ， 可 用 于 挖掘 、 筛 选 和 预测 与 橡胶 生物 合成 相关 的 未 知 基 因 ， 为 研究 橡胶 树 产 量 形 成 的 
分 子 调控 机 理 提供 理论 基础 。 
从 多 个 角度 研究 了 “功能 相关 基因 的 表达 相关 ”这 一 现象 。 在 线虫 (Roy et al, 2002; 
Lercher et al., 2003). ERE (Cohen et al., 2000; 王 晨光 等 ，2006)、 果 蝇 〈Spellman & Rubin, 
2002) 和 人 类 【 郭 政 等 ，2003; Lercher et al., 2002) 等 不 同 物种 上 陆续 发 现 了 染色 体 上 邻近 
基因 表达 相关 这 一 现象 。 酵 母 中 互 作 蛋 白质 的 编码 基因 【王磊 等 ，2004; 欧阳 玉 梅 ，2008; 
王 晨 光 等 ，2006)、 和 蛋白 质 复 合 物 的 编码 基因 (EAF, 2004; Janse et al., 2002) 表达 显著 
> HX, JEHLZUBSEHX XoBR AED GAB ORDERS EE, 2004), EHE A AX E 
LO 于 具有 相同 的 亚 细 胞 定位 〈 欧 阳 玉 梅 ，2008; 王 晨光 等 ,2006) ; 对 基因 数据 库 的 分 析 表 明 ， 
Y 同一 代谢 通路 上 的 基因 表达 倾向 于 高 度 相 关 (Curtis et al., 2005; Harris et al., 2004; Kanehisa 
c et al., 2004) ， 在 酵母 〈 李 传 星 等 ，2004; 王 晨光 等 ，2006) 和 狗 〈 华 琳 等 ，2008) 的 研究 
中 进一步 证 实 了 这 种 相关 性 ;这些 结果 提示 真 核 生物 中 的 基因 存在 模块 化 的 共 表达 趋势 。 
“功能 相关 基因 的 表达 相关 ”有 具 有 实验 条 件 依赖 性 : 不 同 的 实验 条 件 对 这 种 相关 性 有 
c TRI GJG, 2006; Bader et al, 20040 。 橡 胶 生 物 合成 途径 中 也 存在 类 似 的 实验 
Y— 条 件 依赖 性 , 如 在 该 研究 的 S/2D d3 的 割 胶 制度 下 ,MYC7 和 MYC2 显著 正 偏 相关 , 而 在 S/2D 
c2 dl 的 割 胶 制 度 下 则 刚好 相反 (Zhao yue et al., 2011)。 研 究 基因 组 内 和 跨 基因 组 的 基因 表达 
n. 谱 与 蛋白 质 互 作用 之 间 的 相关 性 ， 发 现 与 人 类 、 小 鼠 等 物种 相 比 ， 酵 母 的 蛋白 质 相 互 作用 与 
,一 基因 表达 谱 的 相关 性 最 弱 (Bhardwaj & Lu, 2005); 推测 “功能 相关 基因 的 表达 相关 ”具有 
pS 物种 /进化 程度 依赖 性 。 
“HbCOI1-HbJAZSs-HbMYCSs” 彼 此 间 、 家 族 成 员 间 以 及 他 们 与 FPS、SRPP、REF、HRT2 
"e 等 橡胶 生物 合成 关键 酶 的 相互 作用 ( 马 弗 明 , 2010; 赵 悦 ，2011; XE, 2011; HÆ, 2013; 
FT AZA, 2014; HE, 2015; -E»8, 2016; "f. 2016; 邓 小 敏 , 2018; Deng et al, 2018) X 
明 他 们 在 调控 橡胶 生物 合成 过 程 中 具有 协同 作用 ， 为 这 些 基 因 的 表达 相关 性 提供 了 理化 基 
础 。 


橡胶 树 的 产量 性 状 是 一 种 典型 的 次 生产 量 性 状 ， 由 反复 合理 的 收获 胁迫 《〈 割 胶 ) 形成 ， 
是 一 个 “形成 -收获 -再 形成 -再 收获 ……” 的 过 程 〈 林 位 夫 ，2012) ; 2 次 割 胶 之 间 的 橡胶 再 
生 以 及 相关 基因 的 表达 是 一 个 动态 过 程 ; 由 于 基因 表达 的 动态 性 、 瞬 时 性 、 多 样 性 和 复合 性 ， 
不 同 mRNA 的 稳定 性 和 降解 速率 不 同 ， 可 能 会 影响 特定 实验 条 件 下 的 基因 表达 水 平 
(Kruglyak & Tang, 2000) 。 因 此 ， 系 统 的 追踪 这 些 相 关 基 因 在 排 胶 过 程 中 的 表达 相关 性 是 
今后 的 一 个 研究 方向 。 
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